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Abstract. We report studies on first-order Fermi acceleration in parallel modified shock waves with a large scattering center 
compression ratio expected from turbulence transmission models. Using a Monte Carlo technique we have modeled particle 
acceleration in shocks with a velocity ranging from nonrelativistic to ultrarelativistic and a thickness extending from nearly 
steplike to very wide structures exceeding the particle diffusion length by orders of magnitude. The nonrelativistic diffusion 
approximation is found to be surprisingly accurate in predicting the spectral index of a thick shock with large compression ratio 
even in the cases involving relativistic shock speeds. 
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1. Introduction 

First-order Fermi acceleration in shocks with steplike veloc- 
ity profile is well known to be an efficient way of producing 
nonthermal particle populations with power law (see, e.g., re- 
views of Drury fT955l and Kirk & Duffy fT999l for nonrel- 
ativistic and relativistic studies, respectively). In the nonrela- 
tivistic regime with infinitely thin step-shocks this mechanism 
produces particle momentum distributions f(p) oc p~ s with 
spectral index s depending only on the compression ratio of 
the shock, r, as s — 3r/(r - 1). For the relativistic case also the 
shock speed V\ affects, and the spectral index - in the case of 
isotropic pitch-angle diffusion - can be written as 
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as was recently shown by Keshet & Waxman ( 2005 1. 

In a modified shock - i.e., in a shock with a nontrivial ve- 
locity profile and finite shock thickness - the acceleration ef- 
ficiency drops as the thickness increases (e.g., Schneider & 
KirkQ^Hl Virtanen & Vainio l2"0"0"3l . Again, in the nonrela- 
tivistic case the resulting spectral index as a function of shock 
thickness and compression ratio can have a simple analytical 
solution, found by Drury et al. {I982i, whereas the matter be- 
comes much more complicated in the relativistic regime (e.g., 
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Schneider & Kirk 1989). In their nonrelativistic study Drury et 
al. ([1982) showed that while the produced spectral index tends 
to the well known step-shock limit when the shock thickness 
approaches zero, the high energy part of the spectrum can have 
a power law even when the transition is large compared to free 
path of an accelerating particle. The produced spectral index 
depends on the thickness of the transition region and the com- 
pression ratio as 

3v i /, 1 ^2 \ 3r / 1 \ 

s = — 1 + — = 1 + , (2) 

Vi-V 2 \ /3V l -V 2 ) r-l\ Pir-l)}' 

where f3 is a parameter inversely proportional to the shock 
thickness. For a self-consistent cosmic -ray dominated shock it 
has the value of /3 — ^(1 + y c ), where y c is the adiabatic index 
of the cosmic -ray gas. Later Schneider & Kirk (1989) found 
the nonrelativistic diffusion approximation with spectral index 
of the form of Eq. (0 to agree well even in cases in which the 
thickness of the shock transition is of the order of the parti- 
cle mean free path where the diffusion approximation should 
not be mathematically justified. In the relativistic regime the 
diffusion approximation was, expectedly, found to fail to ap- 
proximate the produced spectral index. One should also note 
that Drury et al. fTW% and Schneider & Kirk ( 1989) neglected 
the second-order Fermi acceleration in their studies. The effect 
of this mechanism in the shock downstream has been studied 
analytically in modified nonrelativistic shocks by Schlickeiser 
(1989) and numerically in relativistic shocks by Virtanen & 
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Fig. 1. Ratio of the scattering center compression ratio to that 
of the gas r for different Alvenic Mach numbers M. The spec- 
tral index of the turbulence power spectrum is q — 2, corre- 
sponding to a momentum-independent particle mean free path. 
M c = yfr. See text for details. 



Vainio (2005 1, and it leads to more complex dependence of 
the produced spectral index of the shock's Peclet and Alfvenic 
Mach number and compression ratio. 

From Eq. it becomes obvious that even shocks with 
large thickness (i.e., small ft) can accelerate particles efficiently 
if only the compression ratio is sufficiently high. This may 
seem irrelevant, as the gas compression ratio of a shock is gen- 
erally limited to rather small values, e.g., r < (y c + l)/(y c - 1) 
for cosmic-ray dominated nonrelativistic shocks. However, the 
compression ratio felt by the accelerated particles is not that of 
the flow; as already noted by Bell ( 1978 1, one needs to consider 
the finite phase speed of the hydromagnetic waves responsible 
for particle scattering around the shock waves. Taking the fluc- 
tuations to be Alfven waves, the problem of turbulence trans- 
mission at the shock is tractable and the effective scattering 
center compression ratio is found to be much larger than the 
gas compression ratio for shocks of low Alfvenic Mach number 
(Vainio & Schlickeiser fT998lfT999l Vainio et al. t20Q3ll20Q5t . 

2. Scattering-center compression ratio 



As shown by Vainio & Schlickeiser ( 1998 1 for nonrelativistic 
and Vainio et al. J2003I I2005l > for relativistic speeds, shocks 
with low to moderate quasi-Newtonian Alvenic Mach number 
M > -\fr transmit the Alfven waves from upstream to down- 
stream so that most of the downstream waves propagate an- 
tiparallel to the direction of the shock-frame plasma flow. Thus, 
the scattering center compression ratio, calculated assuming a 
vanishing upstream cross-helicity, 



Vi(l +//c2Va2V 2 /c 2 ) 
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(3) 



tends to infinity as the Alfvenic Mach number of the shock 
M = u\/u\i (where u\ and mai are the shock proper speed and 
the upstream proper Alfven speed respectively, both measured 
in the rest frame of the upstream plasma) approaches the crit- 
ical Mach number M c = y/r, since V 2 — > Va2 and H c2 — > -1 



at this limit. Here, Va2 and H c2 are the Alfven speed and the 
wave cross helicity in the downstream region, and they can be 
calculated from the shock jump conditions. This is shown in 
Fig.[flfor parallel shocks with three velocities assuming a mo- 
mentum independent scattering mean free path and a vanishing 
upstream cross-helicity. For details on how to compute r^, see 
Vainio et al. (I2003II2005I . 

The scattering hydromagnetic fluctuations have wave- 
lengths comparable to the Larmor radius of the resonant par- 
ticles. Thus, for shock waves discussed in this paper, the 
Alfven waves actually see the shock front as a thick struc- 
ture. Thus, their transmission becomes governed by the WKB 
theory. A detailed analysis of this case is beyond the scope 
of this short Note but qualitatively one expects wave trans- 
mission in these shocks to show very similar behavior as in 
steplike shock waves: antiparallel propagating waves are am- 
plified much more than the parallel propagating ones because 
of the conservation of the wave action, (V ±V\) 2 P ± (f, x)/Va = 
const., where / = (V ± V^k/ln and P +( ~\f,x) are the (con- 
served) wave frequency and the power spectrum of the parallel 
(antiparallel) propagating waves, respectively. For the nonrel- 
ativistic case and for spectral index q of the fluctuations one 
finds after a short calculation 



(4) 



for the downstream intensities of waves propagating parallel 



(antiparallel) to the flow, I^m- For a vanishing upstream cross 
helicity we get 



_ (M + l)i +l (M - V0 9+1 - (M — 1)« +1 (M + v7) 9+I 
c2 ~ (M + iyi +1 (M - yf?Y' +1 +(M - \)i + \M + Vr> 9+1 ' 

Thus, rjfc — > oo as M — > yr, as in our earlier calculations. An 
example of this is shown in Fig. for nonrelativistic shock 
(dash-dot-dashed line, labeled "WKB 0.1c"). Of course, the 
singularity is not reached in practice because the wave pres- 
sure in the downstream region becomes large and makes the 
compression ratio stay below M 2 in a self-consistent calcu- 
lation (Vainio & Schlickeiser 1999 1. Nevertheless, large scat- 
tering compression ratios are expected for low-Mach number 
shocks even for large shock thickness. 



3. Particle acceleration 

We have used test-particle Monte Carlo simulations to calcu- 
late the particle energy distribution N(E) cc E~ a spectral in- 
dices cr — s - 2 resulting from first-order Fermi acceleration at 
the shock for parallel shocks. Detailed description of the code 
is given elsewhere (Appendix A of Virtanen & Vainio 2005 1. 
The simulations follow test-particles in a one dimensional flow. 
Small pitch angle scatterings are performed in the frame co- 
moving with the scattering centers. Instead of using two scat- 
tering wave fields, /*, we calculate the effective wave speed as 
H c Va and apply only a single scattering wave field propagating 
at this speed relative to the plasma. This eliminates stochastic 
acceleration from the model and allows us to study the effect 
of first-order Fermi acceleration at the shock, only. 
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We use the hyperbolic tangent function of Drury et al. 



( 1982 1 for the scattering center speed (denoted hereafter by V) 
profile across the shock: 



V(x) = Vi 



V, -V 2 



1 + tanh 



2c 



A 



(6) 



Here, A is the (constant) particle mean free path and the param- 
eter f3 the same as in Eq. (0. In this study we consider /3 as a 
free parameter and simulate shocks of different widths by vary- 
ing its value. Defining the diffusion length of the particles in the 
usual nonrelativistic manner as di = kAc/Vi, we can write the 
inverse shock thickness as W~ l = \/3{d^ 1 - d 2 l ). Thus, nonrel- 
ativistic simulations with (5 <K 1 correspond to shocks with a 
thickness much larger than the diffusion length of the acceler- 
ated particles. 

For each simulation we fix the proper speed of 
the shock relative to the upstream scattering centers, 
Mi = T\V\ - c(T^ - l) 1 ^ 2 , and the (scattering center) compres- 
sion ratio r/t = V1/V2. Particles are given a small initial energy 
and they are injected into the shock in the upstream so that 
they have time to isotropize before reaching the shock transi- 
tion region. Particles are then traced in the shock until they hit 
the downstream escape boundary at x 2 - This boundary is set 
sufficiently far away from the shock so that the particles have 
enough time to isotropize in the rest frame of the downstream 
scattering centers. A probability of return is calculated and ap- 
plied for each particle crossing the boundary; the method is 
described e.g., by Ellison et al. ( 1990l. Particle splitting is used 
to improve the statistics at high energies. Energy losses due to 
synchrotron emission are not considered in the simulations. 

Simulations were run for shock proper speeds u\ - 
{0.1, 1, 10} c, compression ratios = {3, 10, 30), and for/? hav- 
ing values between ~0.07 and ~ 1000. The resulting spectral 
indices are shown in Fig.|2]together with the predictions of the 
nonrelativistic diffusion approximation from Eq. 0. 

In the nonrelativistic case (Vi ~ u\ — 0.1 c) the results are 
as expected: spectral indices are in accordance with the diffu- 
sion approximation in all cases, even where the shock thick- 
ness is comparable to the diffusion length of the particles. This 
is the case with all compression ratios studied. For the mildly 
relativistic case (u\ = 1.0c, i.e., V\ ~ 0.7c and F\ * 1.4) and 
the lowest compression ratio = 3 (corresponding to the gas 
compression ratio of a relativistic parallel shock) the simulated 
indices start to differ from the diffusion approximation, and 
the differences grow even larger for the ultrarelativistic case 
(mi = 10 c with V\ ~ 0.995 c and Ti * 10). At the step-shock 
limit, where /3 — > 00, the index from the ultrarelativistic shock 
tends to <x ss 2.2 as expected, and also the u\ — c case pro- 
duces indices harder than the diffusion theory predicts. This, of 
course, is expected from previous studies. As the shock thick- 
ness increases (J3 decreases) behavior similar to that reported 
by Schneider & Kirk ( 1989 1 is seen: while for thin shocks the 
diffusion approximation predicts harder spectra than those ob- 
tained from simulation, for shocks with wider transition the 
spectra are softer than the values from the diffusion approxi- 
mation. 

The most interesting result of our study is, however, that 
as the compression ratio increases above the traditional gas- 
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Fig. 2. Particle energy spectral indices as a function of (3 for 
shock proper speeds Mi = 0.1 c (diamonds), 1 .0 c (squares), and 
10.0 c (circles) for scattering center compression ratios = 3, 
10, and 30; the crosses at the right end denote values obtained 
from Eq. {1} for step shocks (J3 — > 00). Solid lines show the 
predictions of the nonrelativistic diffusion approximation. The 
statistical errors of the fitted spectral indices are of the order of 
0.040 or less for all datapoints. 



compression ratio values, the nonrelativistic diffusion approx- 
imation works rather well even for relativistic shocks. Small 
deviations are seen in all cases where relativistic effects are 
present, but the spectral indices can be approximated by Eq. 
to an accuracy of the relevant observations, especially for 
the case = 30. Furthermore, even shocks with very wide 
transition can produce spectra that are hard, compared to those 
usually considered to be produced in parallel shocks. 



4. Conclusions 

Based on our analysis, we arrive at the following conclusions 
of particle acceleration at modified shocks. 

I. Alfven-wave transmission can lead to large scattering 
center compression ratios in parallel shocks regardless of 
the thickness of the shock front if the Alfvenic Mach num- 
ber of the shock is close to the critical value M c = s/r. 
II. Modified parallel shocks can accelerate particles effi- 
ciently even if their thicknesses exceed the diffusion 
length of the accelerated particles, if the scattering cen- 
ter compression ratio is large. 
III. Nonrelativistic theory of particle acceleration in modified 
shocks, Eq. Q, gives a good approximation to the spec- 
tral index at shocks with large scattering center compres- 
sion ratio even for ultrarelativistic shock velocities. 
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